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O ■ ABSTRACT 

■ Combining newly identified and previously known BL Lacs from the RASS-Green 

Bank (RGB) catalog, we present a sample of 127 BLLacs, the largest ever derived from 
^ , a single uniform survey. A Complete sample of 33 objects brighter than 0=18.0 mag 

^ I is also presented. These samples are compared to other known BLLac samples 

and are generally found to exhibit properties intermediate between those of the 
previously disparate classes of High and Low energy peaked BLLacs (HBLs and LBLs, 
respectively). This result is most dramatic in the distribution of the X-ray to radio 
logarithmic flux ratios, where the RGB BLLacs are shown to peak precisely where 
the sharp dichotomy between the two subclasses was previously seen. The a^o vs. Oox 
diagram also shows the RGB sample smoothly bridges the gap between the previously 
distinct subclasses of LBLs and HBLs. The range of broadband Spectral Energy 
Distributions (SEDs) exhibited by the RGB objects also shows that contrary to prior 
claims, searches based on relatively deep surveys cannot limit followup spectroscopy to 
targets with a narrow range of SEDs since BL Lacs clearly constitute a homogeneous 
population with a wide range of SEDs. 

Similar to results based on the EMSS and 1 Jy BLLac samples, we find a weak, but 
statistically significant correlation between the composite X-ray spectral index Oxox 
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and the radio-optical spectral index a^o- This implies that the more LBL-like RGB 

BL Lacs have a secondary source of X-ray emission, possibly from an inverse Compton 
component. This result, in addition to other characteristics of the RGB sample, 
indicates that the simple unified scheme which postulates HBLs and LBLs differ solely 
by orientation may be in need of revision. We also present both the X-ray and radio 
logN— logS distributions for which the competing HBL/LBL unification scenarios have 
differing predictions. The unknown effects of the triple fiux limit inherent in the RGB 
Complete sample makes quantitative analysis uncertain, but the characteristics of the 
RGB sample compare well both with results obtained from previous samples and with 
general theoretical predictions based on a simple Monte Carlo simulation. Our analysis 
indicates that the unimodal distribution of BL Lac properties found in the RGB sample 
likely reliably reflect the underlying population, while the bimodal distribution found 
in earlier studies arose primarily from observational selection effects. The presence of 
not only intermediate, but also extreme HBL and LBL objects is the RGB survey's 
unique strength and offers clear avenues for future studies which can undoubtedly 
address the question of how HBLs and LBLs are related. 

Subject headings: BL Lacertae objects: general — galaxies: active — radio continuum: 
galaxies — surveys — X-ray: general 
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1. Introduction 

BLLacs comprise a rare subclass of Active Galactic Nuclei (AGN) and are characterized 
by a lack of prominent emission lines, a highly variable nonthermal continuum and strong, 
variable optical polarization (KoUgaard 1994; Urry and Padovani 1995). Additional characteristics 

include the lack of a UV-excess (or "blue bump") and a core-dominated radio morphology 
(Angel and Stockman 1980; Wardle et al. 1984; Laurcnt-Muehleisen et al. 1993). The dominance 
of nonthermal radiation at all wavelengths ranging from the radio to gamma ray regimes is 
well-established (Sambruna et al. 1996, and references therein) and makes BL Lacs particularly 
interesting laboratories for the study of AGN phenomena. 

Numerous studies have shown that BLLacs contain relativistic jets which produce narrow 
cones of beamed emission which makes the observed radiation sensitive to orientation of the jet 
axis relative to the line-of-sight (see Kollgaard 1994; Urry and Padovani 1995, and references 
therein). BLLacs are associated with those objects that are oriented such that their jets lie close 
to the line-of-sight while the parent population of misaligned objects are postulated to be low 
luminosity Fan ar off- Riley Type I radio galaxies (Fanaroff and Riley 1974; Browne 1983). This 
association between BL Lacs and FR I galaxies is one of the stronger links in the "unified scheme" 
of AGN where observed properties are primarily a result of orientation rather than intrinsic 
astrophysical differences (e.g., Antonucci 1993). 

However, it has not been clear whether all BLLac properties can be attributed to 
orientation differences, nor whether BL Lacs constitute a homogeneous class. For example, the 
broadband spectral energy distributions (SEDs) of BL Lacs discovered in X-ray and radio surveys 
differ significantly which has led to the sub classification of BLLacs into X-ray— selected and 
radio— selected objects (XBLs and RBLs, respectively). This has recently been supplanted by 
a new classification "High energy peaked BL Lacs" (HBLs) and "Low energy peaked BL Lacs" 
(LBLs) based on the ratio of X-ray to radio flux densities, Sx/Sr (Giommi and Padovani 1994; 
Padovani and Giommi 1995a). Generally, XBLs tend to be HBLs and exhibit less extreme 
properties than RBLs which are usually LBLs. 

The two BL Lac subclasses exhibit systematically distinct properties, including the degree of 
radio core dominance, optical polarization fraction and duty cycle, fraction of optical host galaxy 
light, and perhaps even parsec-scale jet speeds, megaparsec-scale clustering properties and host 
galaxy optical and radio luminosities (Laurent-Muehleisen et al. 1993; Perlman and Stocke 1993; 
Jannuzi et al. 1994; Wurtz et al. 1996; Kollgaard et al. 1996a, 1999). Many of these characteristics 
are consistent with the unified scheme paradigm if HBLs are objects which lie further from the 
line-of-sight than LBLs. However, orientation by itself cannot explain the apparent lack of BL Lacs 
with properties intermediate between the LBL and HBL subclasses. 

There is also increasing concern that the simplest unified scheme may not account for all 
subclass distinctions. Intrinsic as well as orientation differences may be present. For example, 
estimates of jet speeds and angles to the line-of-sight of HBLs, LBLs and FRI radio galaxies 
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do not appear to be able to account for all differences in the SEDs (Sambruna et al. 1996). In 

addition, estimates of space densities are inconsistent with orientation values, and the cosmic 
evolution of the classes appears to be incompatible (Rector et al. 1999a). 

These issues have proved difficult to address because existing BL Lac samples are still relatively 
small and were generated from shallow surveys which contain only the very brightest objects in 
either the radio or X-ray wavebands. These and other selection effects have produced samples 
biased towards the most extreme HBLs or LBLs with few transitional objects. The sample of 
BL Lacs presented here was generated from a cross-correlation of a deep radio (Laurent- Muehleisen 
et al. 1997) and X-ray catalog and contains BL Lacs with the full range of properties from HBLs 
to LBLs. This RASS-Green Bank (RGB) BL Lac sample is the largest BL Lac sample yet created 
from a uniformly defined set of criteria. It consists of 127 objects drawn from a correlation 
of the ROSAT All-Sky Survey (RASS) and a reanalysis of the 1987 Green Bank 6 cm radio 
survey (GB96, Gregory et al. 1996). The design and followup spectroscopic observations of this 
sample are presented in Laurent-Muehleisen et al. (1998, hereafter Paper I). Here we concentrate 
specifically on the RGB BLLacs and on what they reveal about the relationship between the 
BL Lac subclasses. 

This paper is organized as follows. In §2 we briefly review the RASS/GB correlation and our 
followup VLA and optical observations. The RGB and "RGB Complete" samples are presented 
in §3 and §4. Sections 5 and 6 analyze these samples' bulk characteristics and discuss their 
astrophysical implications. We assume throughout Ho=100 kms~^Mpc~^, qo=0.5 and define 
spectral indices, a, such that Si, oc 

2. Selection of Candidate Objects 

Both radio and X-ray surveys have proven to be a rich source of new BLLacs. The 
largest purely radio— selected sample is that based on the 1 Jy survey while the most prominent 
X-ray— selected samples are those based on the Einstein Extended Medium Sensitivity Survey 
(EMSS; Gioia et al. 1990; Stocke et al. 1991; Maccacaro et al. 1994) and the HEAO-1 Large 
Area Sky Survey (Wood et al. 1984; Kollgaard et al. 1996b). As the flux limit and/or sky 
coverage of radio and X-ray surveys has improved and the ability to fully identify all objects 
in these surveys has become impractical, the technique of selecting candidate BL Lacs based 
on the broadband SEDs of previously known BLLacs has proven to be highly efficient. The 
Einstein Slew Survey (Elvis et al. 1992; Schachter et al. 1993a; Perlman et al. 1996a), Hamburg 
Quasar Survey (Hagen et al. 1995; Nass et al. 1996), Deep X-ray Radio Blazar Survey samples 
(Perlman et al. 1998) and the optical polarization sample of Kock et al. (1996) were all created by 
spectroscopically classifying sources with counterparts detected concurrently in the more than one 
band. Nevertheless, the number of objects in each of these samples has remained relatively small 
and their selection effects are often difficult to assess (see, e.g., Browne and Marcha 1993). The 
need for a large sample with a minimal number of simple selection criteria is clear. 
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The potential of the ROSAT All-Sky Survey for creating just such a BL Lac sample has been 
noted (Stocke et al. 1989) and followup programs are confirming this prediction. Some programs 
based on the RASS have used the optical polarization or narrowly defined broadband SEDs to 
select BLLac candidates (Kock et al. 1996; Nass et al. 1996). These candidates were chosen 
independent of any knowledge of the radio flux of the source, but followup observations have 
shown that all are detected at centimeter wavelengths (Nass et al. 1996), a result in agreement 
with the assertion that radio-silent BLLacs are very rare or possibly nonexistent (Stocke ct al. 
1990). Our RGB BLLac sample was therefore constructed using previously detected radio 
emission as a criterion. This sample is therefore triply flux-limited (radio, optical and X-ray) but 
imposes no other selection criteria other than location in the northern hemisphere {0° <5<75° ; 
Laurent-Muehleisen et al. 1997). This method coupled with high sensitivity of the RGB survey in 
these three wavebands detects BL Lacs with a variety of broadband SEDs. 

The initial RGB catalog consisted of sources whose positions differed by less than 100" 
in the RASS and a point source catalog created from the 1987 Green Bank 6 cm radio survey 
(GB96, Gregory et al. 1996; Neumann et al. 1994). This new GB catalog consists of 3cr and 
greater confidence sources and has a flux density limit of ~15 mJy in the declination range 
from 30° — 75° and increases to ~24mJy at low declinations (Neumann et al. 1994). In order 
to eliminate spurious RASS-GB coincidences and to determine positions accurate enough for 
reliable optical identification, the 2,127 sources in the RASS-GB correlation were observed at 
high resolution with the NRAO's^ VLA. Two radio catalogs were produced, the first consists of 
1,861 sources for which subarcsecond positions and core radio flux densities were obtained; the 
second consists of 436 sources for which only low resolution data (~8" positional accuracy) were 
obtained (Laurent-Muehleisen et al. 1997). Our analysis showed that all sources whose radio 
and X-ray positions differed by less than 40" are true matches to a high degree of confidence. 
Additional details on the RASS-GB correlation and the followup radio observations can be found 
in Laurent-Muehleisen et al. (1997) and Laurent-Muehleisen (1996)^. 

Optical counterparts were determined via Automatic Plate Measuring scans of the high 
Galactic latitude (>25°) POSSI photographic plates (Kibblewhite et al. 1984). Optical 
counterparts within 3" of RGB sources were identified and both the O (blue) and E (red) 
magnitudes measured (Laurent-Muehleisen 1996; Brinkmann et al. 1997). A looser criterion of 5" 
was used for sources in the low resolution VLA catalog. Spectra were obtained for 169 optically 
bright (0<18.5mag) objects which lacked spectroscopic classifications (Paper I). 

Table 1 summarizes the steps involved in the creation of the final RGB catalog. A discussion 
of the broadband multiwavelength properties of the entire RGB sample can be found in Paper I 
and Brinkmann et al. (1997) as well as Laurent-Muehleisen (1996). 



^NRAO is operated by Associated Universities, Inc., under cooperative agreement with the National Science 
Foundation. 

^This Ph.D. thesis is available via the WWW at http://www-igpp.llnl.gov/people/ slauren.html 
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TABLE 1 
RGB Catalog Creation Summary 

Step # of Sources 



Initial RASS-GB Correlation 2,127 

Followup VLA Radio Sources 2,297 

Reliable RASS-VLA Matches 1,567 

Reliable RASS-VLA-POSS I Matches 971 

Spectroscopically Classified Sources 594 

RGB BL Lacs 127 



Notes 

Match Criterion: Ar/x = 100" 
Some GB sources became >1 VLA source 
1,861 from high radio resolution data 
436 from low radio resolution data 
83 Empty fields 
2 Not observed 
Match Criterion: Avla/x = 40" 
796 from high radio resolution data 
Match Criterion: AvLA/Opt = 3" 
175 from low radio resolution data 
Match Criterion: AyLA/Opt = 5" 
Includes 140 new IDs from Laurent-Muehleisen et al. 1998 
548 from high radio resolution data 
46 from low radio resolution data 
100 are unambiguous BL Lacs (27 are transitional objects) 
107 IDs are from high radio resolution data 
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We noted in Paper I that the operational definition of a BLLac has changed since that 
of Stein et al. (1976) who defined a BLLac as an AGN having a highly variable, linearly 
polarized, nonthermal continuum without optical emission lines. More recently, Stockc ct al. 
(1991) defined a BLLac as an AGN with emission lines whose equivalent width does not exceed 
5 A and whose Ca II H&K breaJi strength (Br4ooo) is <25%^. However, recent observations, 
particularly of the less extreme HBLs, has shown that even this definition is too stringent and 
excludes objects which otherwise exhibit BL Lac-like properties and should therefore be classified 
as such (Marcha et al. 1996; Scarpa and Falomo 1997). Establishing a set of criteria which define 
the BLLac class is particularly difficult because much of the observed emission is contributed 
by orientation-dependent beamed radiation while the creation of unbiased samples requires that 
the classification of objects be based on intrinsic (not observed) characteristics. In Paper I, we 
adopted the classification criteria proposed by Marcha et al. (1996) to distinguish BLLacs from 
ordinary quasars, Seyfert, radio and elliptical galaxies. Briefly, our classification scheme is: 

• If the spectrum is featureless or the only features observed are emission lines with Wa<5A 
(rest frame), the object is classified as a BLLac. 

• If absorption features are present and Br4ooo<25%, we classify the object as a BLLac, 
provided any emission lines present have W;^<5A. 

• If the Ca II break contrast is between 25%— 40%, we classify the object as a possible BLLac 
if any emission line present also has an equivalent width smaller than that required by the 
Marcha et al. (1996) criterion for that particular break strength (see Paper I). 

• If the Ca II break contrast was >40% we classified the object as a galaxy because of the lack 
of spectroscopic evidence for an AGN, although we find no clear discontinuity in properties 
which distinguishes BLLacs from galaxies (Paper I, Figure 3). 

Our spectroscopic observations combined with these criteria produced a sample of 53 RGB 
BLLacs, 38 of which were newly discovered. We now combine these objects with previously 
identified BL Lacs in the RGB catalog and discuss the full sample. 

3. The RGB Sample of BL Lacs 

The RGB sample consists of 127 sources of which 100 are definitive BL Lacs. Many of the 
27 objects which are only probable BLLacs have break contrasts larger than 25%, but adhere to 
the Marcha Br4ooo— criteria. Twenty of the RGB BLLacs belong to the low radio resolution 



^The Ca II break contrast refers to the relative depression of the continuum blueward of the Ca II H&K hues 
(3933 A & 3968 A). 
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subset of the RGB catalog. Because core radio flux characterizes many of the beaming properties 

of BLLacs, we exclude these "low radio resolution objects" from all further analysis, reducing the 
sample of objects we discuss to 107 RGB BLLacs which is the largest sample of BLLacs ever 
cataloged from one survey. 

For completeness, all 127 sources are presented in Table 2. The columns give the (1) RGB 
Name; (2) alternate Name; (3) and (4) J2000 coordinates; (5) RGB 5 GHz core radio flux density 
(in mJy); (6) X-ray flux (in 10"^^ ergs'^cm-^ in the 0.1-2.4keV ROSAT band); (7) O magnitude 
from the APM POSS-I catalog; (8) ROSAT PSPC X-ray spectral energy index; (9) and (10) the 
ttro and Oox spectral indices; (11) redshift and (12) references. As discussed in Brinkmann et al. 
(1997), Laurent-Muehleisen et al. (1997) and Laurent-Muehleisen et al. (1998), the uncertainties 
on the multiband fluxes can be taken to be approximately 20%, 0.5 mag and 25% in the radio, 
optical and X-ray, respectively. 

Flux densities in Table 2 have not been K-corrected, but the two-point spectral indices 
(cKro and aox) are corrected assuming a flat (ar=0.0) radio spectral index and an X-ray spectral 
index as given in column 8, or the mean X-ray spectral index of ax=1.2, valid for a large sample 
of RASS-detected BLLacs (Brinkmann et al. 1997). When a measured redshift is unavailable, 
the K-correction is based on the median RGB BLLac redshift, z=0.16. Although this probably 
underestimates these objects' true redshift (objects with z<0.16 would likely show spectral 
evidence for the host galaxy and hence would have a measured redshift), we feel it is important 
to apply some type of K-correction to all objects since the K-correction always produces a net 
flattening of the aro and ctox spectral indices (given our assumptions). However, the effect is small, 
typically producing an error of <0.05 in either index unless the true redshift is >1.0. 

The optical magnitudes are equivalent O band values, converted from other bands assuming 
«opt=l-0 when no O or B magnitudes were available in the literature. (The notes to Table 2 
also give the original magnitude and band.) If any galaxy absorption features were detected, we 

corrected the magnitudes to include only the contribution from the AGN component using the 
method described in Paper I. Wc note that these corrections have only been made for the newly 
identified objects presented in Paper I because the correction requires the determination of the 
Ca II break contrast which is usually not reported in the literature. 

The K-correctcd radio— optical and optical— X-ray spectral indices arc defined to be 

(gcorc \ / S \ 

-^-^j and aox=0.351 log (^-^ ) and are given in columns 9 & 10. The 

monochromatic ROSAT X-ray flux density was converted from the fluxes given in Table 2 

according to: 



Sx = 4.138 X 10-^^ ■ Fx ■ E 



1 - 



_2.4(i-ax) _ 0.l{i-ax) 

where E is set at 2keV, Fx is in ergs~^cm~^ and Sx is in ergs~^ cm~^keV~^. The values of 



(1) 



-9- 



CKro and Q!ox are accurate to ±0.1 and ±0.2, respectively, taking into account the observational 

uncertainties in the fluxes. This docs not, however, include the effects of real temporal variability 
in these nonsimultaneous multiband data. Generally speaking, variability is unlikely to change the 
flux values by more than 25% and the two-point spectral indices given in the table are therefore 
reasonably representative. 

Finally we note that a few RGB sources identified as BL Lacs in the literature are missing 
from Table 2 and therefore technically not members of the RGB BL Lac sample. RGB J1058+564, 
RGB JlllO+715, and RGB J1610+671B have spectra that were presented in Paper I, but violate 
the radio/X-ray coincidence criterion of Arx<40". The previously known BLLacs RGB J0738+177 
(PKS 0735+178) and RGB J1508+271 are respectively 3.1" and 3.8" from bright optical sources 
and therefore do not appear in Table 2 since they violate our radio/optical coincidence criterion. 
However, both are likely true RGB BLLacs^. 

4. The RGB Complete Sample 

The RGB sample was constructed without imposing any selection criteria other than the 
presence of an optical counterpart within 3" (5" for the low radio resolution sample) and a RASS 
source within 40" of a GB96 radio source. These well-defined criteria allowed us to create the 
"RGB Complete Sample" which consists of optically bright objects (total O magnitude <18.0) 
observed over 3970 deg^ of the sky with a completeness of 94% (Figure 1). 

Our complete survey area contains a total of 183 RGB sources, most of which are emission-line 
AGN (68%) and galaxies (12%). But 33 are optically bright BLLacs (including three candidate 
objects with Call break contrasts between 29-39%) which constitute the RGB Complete BLLac 
Sample. These objects are listed separately in Table 2. The RGB Complete sample is therefore 
flux-limited in three bands: radio, optical and X-ray. The constraints imposed by the RASS survey 
(whose flux limit varies with ecliptic latitude and iVn) and the GB survey (whose flux density limit 
varies slightly with declination) both aflFect the catalog's overall completeness but in a well-defined 
manner. The effects of the optical flux limit are more serious than those imposed by the X-ray 
and radio limits. Assuming all BLLacs have color-color indices in the range 0.1<aro<0.8 and 
0.5< Q!ox<2.1 (see Figure 3), a limiting optical magnitude of O — 18.0 mag implies that the RGB 
Complete sample is only truly complete above Sr=3 Jy and Fx=10~-'^^ergs~^cm~^, criteria which 
are satisfled by none of the BL Lacs in this sample. However, the situation is less dire in practice 
since the RGB Complete sample is constructed from surveys with well-deflned flux limits but 
without any additional selection criteria applied to it. Any additional constraints (such as limiting 



^Source 0930-1-4950 is a BLLac but was incorrectly reported in Laurent-Muehleisen et al. (1997) as belonging to 
the RGB catalog of X-ray and radio-emitting AGN. We therefore also exclude this object from the current paper. 
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Fig. 1. — The 3970 deg from which the RGB Complete sample was selected. The Galactic plane 
(b<25°) and northern extent of the GB survey ((5=75°) are labeled. The region designated by 
the horizontal lines constitutes the RGB Complete Survey area which consists of the following six 

regions: (1) 6'^30"'< a <10''31'^ and 57°40'< 5 <75° (2) 10'^31™< a <12'^48™ and 50°< 5 <75° 
(3) 12'^48'>^< a <15'^ and 40°30'< 6 <75° (4) 15'^< a <19'^ and 15°< 5 <75° and (5) 
2l'^30°'< a <23^30™ and 3°< 6 <34° with the added constraint that b>25°. 

candidates to objects with particular SEDs), however, would be much less useful for characterizing 
the full BL Lac population. 

We note that three previously known BL Lacs nearly, but do not exactly, satisfy the criteria 
for inclusion in the RGB Complete sample. First, lES 2326+174 is 16.8 mag Slew Survey BLLac 
with 27mJy radio flux and 1.5 x 10~^^ erg cm~^ X-ray flux (Perlman et al. 1996a). This 
source lies in a region of the All-Sky Survey where the standard data screening software flnds 
very little acceptable data, reducing the effective exposure time to essentially zero. This source 
therefore did not appear in our original RASS— GB correlation and is therefore missing from the 
current RGB BLLac sample. Second, as mentioned in §3, RGB J1058+564 and RGB 1508+271 
are BLLacs that respectively violate our radio/X-ray and radio/optical offset criteria. With these 
possible exceptions, our 33 source optically bright RGB BLLac sample is complete given our 
selection criteria and therefore constitutes a useful sample for statistical study. 
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5. Properties of the RGB Sample 

To assess the characteristics of the RGB sample, we compare it with four large, well-defined 
samples of BLLacs currently available in the literature: the 1 Jy, HEAO, EMSS and Einstein Slew 
Survey samples (Stickel et al. 1991; Stocke et al. 1991; Perlman et al. 1996a; Remillard et al. 1986; 
Schwartz ct al. 1989). We will treat the RBL 1 .Jy sample as representing the range of properties 
associated with LBLs (Low energy peaked BLLacs), although the recent assertion that some 
IJy objects may be misclassified microlensed sources should be noted (Stocke and Rector 1997). 
The XBL EMSS and HEAO samples will be combined to represent HBLs (High energy peaked 
BLLacs). As will be shown below, the Slew Survey objects tend to exhibit properties intermediate 
between HBLs and LBLs and are in that sense similar to the RGB sample. The occassional 
duplication of a source in more than one sample is ignored here. 

We compare the distributions of redshift, three representations of SED shape — Sx/Sr ratio, 
Q!xox=Q!ox— cnx (Sambruna et al. 1996) and location in the a^o vs. Oox plane — and the radio and 
X-ray BLLac logN— LogS relations. Data for the comparison samples were obtained from the 
literature, mainly from Laurent-Muehleiscn et al. (1993), Kollgaard ct al. (1996b), Sambruna 
et al. (1996), Perlman et al. (1996a), Perlman et al. (1996b), Urry ct al. (1996), Brinkmann et al. 
(1997) and Paper L The median values for the various samples are given below in Table 3 and 
were calculated using the Astronomy SURVival analysis software (ASURV, Rev. 1.2, LaValley 
ct al. 1992, which can be obtained from littp://www. astro. psu.edu/statcodes) which properly 
handles the upper and lower limits present in the data. (One source has an upper limit to its core 
radio flux density, RGB J1000-I-225A.) Median values are calculated using the maximum-likelihood 
Kaplan-Meier estimator and the statistical significance of any differences between two samples is 
estimated with the logrank and Peto & Peto generalized Wilcoxon tests (Feigelson and Nelson 
1985). 



Table 3. Median Properties of BL Lac Samples 



Property 



HEAO EMSS 



XBL 

(HEAO-F 
EMSS) 



Slew 



RGB 
Complete 



RGB 1 Jy 



Redshift 0.12 0.30 0.20 
log Sx/Sr -4.78 -4.84 -4.79 
axnx -0.25 -0.13 -0.17 



0.16 0.13 
-4.56 -5.56 
-0.23 -0.11 



0.16 0.50 
-5.61 -6.89 
-0.14 0.23 
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Fig. 2. — The redshift distribution for the various BL Lac samples, a. The XBL sample of BL Lacs 
which consists of both the EMSS and HEAO objects. The hatched region denotes the HEAO 
objects and the top solid line is the sum of both samples, b. The Slew survey sample of BL Lacs, 
c. The RGB sample of BL Lacs. The hatched region denotes the RGB Complete sample of bright 
objects, d. The IJy sample of RBLs. The x's on the far right represent objects with unknown 
redshifts. 
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5.1. Redshift Distribution 

Figure 2 shows the distribution of redshifts for the RGB sample as weU as the XBL, Slew 
Survey, and 1 Jy samples. A typical HBL clearly resides at a lower redshift than a typical LBL, 
a trend which is predicted by the unified model of Fossati et al. (1997) which asserts that HBLs 
constitute the intrinsically lower luminosity sources. Although redshifts are known for only 
approximately half (59) of the RGB sources, RGB redshifts span nearly the entire range exhibited 
both by HBLs and LBLs but they are heavily weighted toward lower redshifts with a smaller 
median redshift (0.16) than most samples. The median redshift of the RGB Gomplete sample is 
essentially the same as the full RGB sample. 

Not surprisingly, many of the lowest redshift objects in the RGB BL Lac sample exhibit a 
Ca II break contrast: 13 of the 17 RGB BLLacs with measurable Ca II break contrasts lie below 
the median redshift. This implies that earlier samples which required Ca II break contrasts to be 
less than 25% are missing a significant fraction of low redshift objects, a result also obtained in the 
preliminary work on the "REX" BLLac sample (Caccianiga et al. 1999). Perlman et al. (1996a) 
noted the typical redshift of a Slew Survey BLLac is much lower than that of an EMSS object, 
perhaps indicating the EMSS sample is incomplete due to misclassification. Recent observations 
indeed show that four EMSS objects should be reclassified as BLLacs based on the Marcha 
et al. (1996) criteria (Rector et al. 1999a). These points illustrate the ambiguity in the current 
definitions for BLLacs. The Marcha et al. (1996) criteria, adopted in our RGB sample, are 
however an important step toward a standard, meaningful definition that is largely independent 
of orientation and encompasses the lower luminosity (and presumedly more numerous) objects. 

5.2. The aro vs. aox Diagram 

In Figure 3 we present the a^o vs. a^^ color-color diagram for the RGB and comparison 
BL Lac samples. The spectral indices for the comparison samples have been recalculated using our 
assumptions (§3) and based on data in the literature. When the redshift is unknown, we use the 
median redshift for objects in that sample (see Table 3). However, we were not able to correct the 
optical fluxes from other samples to include emission from only the AGN component, as we did 
for the RGB sample because the Ca II break contrasts are generally unavailable. The correction 
for host galaxy starlight is typically an increase (decrease) in a^o (o^ox) of 0.06 (0.11) for the 17 
"corrected" RGB objects which has little effect on the Oro— "ox diagram. Additionally, unlike most 
diagrams of this type, the radio flux density used to calculate the aro spectral index includes the 
5 GHz flux from only the arcsecond-scale radio core which best represents the beamed component. 
This is an important point, as extended emission is frequently comparable to core emission in 
HBLs (Laurent-Muehleisen et al. 1993; Perlman and Stocke 1993). 

Figure 3 shows the RGB BLLacs exhibit a smooth distribution in Oox; ranging from 0.46—1.72, 
with no hint of the bimodality which has been previously widely discussed (e.g., Brinkmann et al. 
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Fig. 3. — The a^o vs. aox color-color diagram for various samples of BL Lacs. The flux densities 
used to calculate the spectral indices have all been K-corrected and converted to 5 GHz (radio core) , 
4400 A, and 2 keV for the three bands. Objects common to more than one sample are plotted using 
the symbols representing all samples to which they belong. The EMSS-defined class boundaries 
are also shown for comparison. Note that the RGB sample spans both the traditional HBL region 
(defined by the EMSS and HEAO samples) and the LBL region (defined by the 1 Jy sample). 

1996). This result is in agreement with early results from the "DXRBS" (Perlman et al. 1998) and 
"REX" (Caccianiga et al. 1999) samples which show that a large number of intermediate objects 
exists and that there is no clear separation of the HBL and LBL subclasses (see also §5.3 below). 
This implies the apparent bimodality could have been caused by selection effects inherent in the 
two best-studied, previously known samples: the EMSS and 1 Jy, a question we will examine in 
greater detail in §6. 

Both the RGB and Slew survey objects lie in regions of flatter a^-o and steeper aox than 
the EMSS or IJy BLLacs. Generally, the RGB BLLacs lie along a horizontal band defined by 
0.2< ckro<0-6. All but one of the RGB objects with Q;ro>0.6, also have aox>l-0 and therefore 
have LBL-like SEDs, although few RGB objects appear to be as extreme as the majority of 1 Jy 
BLLacs. This likely occurs because any RGB object fainter than 0!^18.5mag (the magnitude 
limit for most of the RGB sample) with Q!ro>0.6 must be very radio bright, Sr;^200mJy, a flux 
density large enough to exclude the majority of RGB objects. In contrast, the 1 Jy sample has an 
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optical magnitude limit of V~20mag and therefore contains objects with very steep values of CKro- 

As the identification of the RGB catalog (Laurcnt-Muehleisen et al. 1997) is extended to fainter 
optical magnitudes, we expect objects with flat a^^ but steep Oro spectral indices (i.e., extreme 
LBL-like BL Lacs) will be discovered. 

5.3. X-ray to Radio Flux Density Ratios 

Previous BL Lac samples show a clear bimodality in the ratio of the X-ray to radio flux 
densities of HBLs and LBLs at log Sx/SrC±— 5.5 (Padovani and Giommi 1995a; Perlman et al. 
1996a; Brinkmann et al. 1996). In contrast, the RGB BLLacs have a median X-ray to radio 
logarithmic flux density ratio of —5.61 (Table 3). The distribution of the flux ratios (Figure 
4) shows no evidence for a sharp division between the two subclasses. The RGB catalog is 
therefore the first to contain large numbers intermediate BLLacs, although hints that these 
objects existed have been previously reported by Perlman et al. (1998), Caccianiga et al. (1999) 
and Laurent-Muehleisen et al. (1998). 

Statistical two-sample tests show that the RGB and the RGB Complete samples are consistent 
with having been drawn from the same distribution. This indicates that the bright optical 
magnitude limit imposed on the RGB Complete sample does not drastically affect the overall 
distribution of Sx/Sr and, by extension, the fraction of HBLs vs. LBLs, a property for which the 
different unification models have significantly different predictions. Studies of the properties of 
the Complete Sample can therefore yield important insights into the origin of the HBL vs. LBL 
subclasses. 

5.4. High Energy Continuum: axox 

The X-ray emission mechanism in both HBLs and LBLs is most likely dominated by 
synchrotron radiation at low energies and inverse-Compton (IC) processes at higher (7-ray) 
energies (Blandford and Rees 1978; Ghisellini et al. 1985; Bregman et al. 1987). The shape of 
the X-ray spectrum is useful for determining at what energy this transition takes place which is 
important for understanding the overall energy budget and underlying jet physics in BL Lacs. 

While X-ray spectral indices characterize the high energy continuum, the composite 
X-ray/optical spectral index, axox=CKox — ckx more precisely measures changes in the SED between 
the optical and soft X-ray bands by distinguishing the relative importance of IC and synchrotron 
emission processes. If q;xox<0, then the X-rays lie along a powerlaw or steepening synchrotron 
continuum. A positive value of axox represents a concave spectrum and is likely caused by a hard 
IC component. 
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Fig. 4. — The distribution of the logarithm of the X-ray to radio flux density ratios for various 
samples of BL Lacs a. The XBL sample b. The Slew survey sample c. The RGB sample d. The 
1 Jy sample. Hatched regions are as described in Figure 2. Unlike many other samples, we consider 
here only the core radio flux density. The division between HBLs and LBLs, as defined by the flux 
ratios of previously known samples, clearly occurs at log Sx/Sr=— 5.5. The RGB sample shows no 
such dichotomy and also includes objects with traditional HBL- and LBL-like SEDs. 
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Fig. 5. — The distribution of the difference between the optical X-ray and X-ray spectral indices, 
Q^xox=«ox— cKx5 for various BLLac samples a. The XBL sample b. The Slew survey sample c. The 
RGB sample d. The 1 Jy sample. Hatched regions are as described in Figure 2. 
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Fig. 6. — The Oj-o vs. axox diagram for the RGB BLLacs. Fihed circles represent LBLs and X's 
represent HBLs. There is a weak (P=92%) correlation present. Note that LBLs tend to be the 
objects with the steepest a^o and positive ("convex") axox = ckox — ckx. 



If LBLs and HBLs differ solely by orientation, then it should be possible to explain differences 
in their SEDs by invoking beaming models and varying only one free parameter, the angle to 
the line-of-sight. However, Sambruna et al. (1996) find that there are differences in the SEDs of 
LBLs and HBLs which are not attributable to orientation alone: objects with steeper a^o spectral 
indices (generally LBLs) have a tendency to also have more positive values of a^ox- This indicates 
the presence of a second, presumably hard IC, X-ray component and suggests there are intrinsic 
differences between LBLs and HBLs which are independent of orientation. 

Figure 5 shows the distribution of a^ox for the XBL, Slew, RGB and 1 Jy samples. Only 
those objects with measured X-ray spectral indices are included (82% of the RGB sample). As 
seen for other parameters, the RGB sample spans the range exhibited by the extremes of LBLs 
and HBLs and many sources exhibit intermediate properties; the RGB Complete sample exhibits 
similar characteristics (Table 3). To test the results of Sambruna et al. (1996), we searched for a 
correlation between a^-o and axox in the RGB sample (Figure 6). Objects have been divided into 
HBL- and LBL-like classes based on their X-ray to radio flux ratios. There is a correlation with 
moderate statistical significance (P~92% using a nonparametric Spearman's p statistic) for the 
RGB BLLacs. This is weaker than the correlation found by Sambruna et al. (1996) who found 
P>99.99% using the EMSS and 1 Jy samples. However, if we supplement the intermediate RGB 
objects with BLLacs from the EMSS, HEAO and 1 Jy samples, the probability that a correlation 
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is present increases to >99.99%. It is not surprising that the addition of the extreme HBLs (which 
He in the bottom left corner of the diagram) and the extreme LBLs (which He in the upper right 
corner) strengthen the correlation. This shows that the full range of high energy spectral shapes 
exhibited by BL Lacs does correlate with BL Lac subclass, confirming the results of Sambruna 
et al. (1996). 

5.5. X-ray logN— logS Relationship 

Both the X-ray and radio logN-logS distributions of the EMSS and 1 Jy BL Lacs have been 
studied extensively. However, as the results presented here show that BLLacs do not belong 
to two distinct subclasses, conclusions drawn from studies of the logN-logS distribution of only 
the two extremes of the BLLac population may be misleading. Nevertheless, theoretical models 
based on the 1 Jy and EMSS data have made interesting predictions. Based on the standard 
orientation-based beamed jet model unifying FRI radio galaxies and BLLacs, the known FRI 
luminosity function and observational constraints on jet Lorentz factors, Urry et al. (1991) 
calculate the radio radio logN-logS distribution of BL Lacs down to 1 mJy. We compare the 
Complete RGB sample to these predictions and other BL Lac samples. 

One caveat regarding the completeness of the sample bears mentioning. We characterize the 
RGB BL Lacs' radio emission using core flux densities which can often fall significantly below the 
total flux density as measured in the original Green Bank survey. As a consequence, only those 
objects with low radio core-to-lobe ratios populate the lowest radio flux density bins. (Faint radio 
sources with high radio core-to-lobe ratios will have total radio flux densities below the GB survey 
flux limit and will be missing from our survey.) Only the lowest radio flux density bins (<20mJy 
based on the GB96 survey flux limit) are affected by this source of incompleteness. 

Figure 7 shows the X-ray logN-logS distribution for the RGB Complete, EMSS and Hamburg 
Quasar Survey samples (HQS, Hagen et al. 1995; Nass et al. 1996). Here we truncate the EMSS 
and HQS samples so that they match the selection criteria present in the RGB Complete catalog, 
namely the objects must be brighter than 0=18.0 mag in the optical and bright enough in the 
radio to be present in the GB96 catalog. This eliminates 19 of the 22 EMSS objects, while the 
HQS sample is diminished from 61 objects to 27. 

Also shown in the figure is the Euclidean no-evolution model defined by N(>Fx)ocFj^' and 
arbitrarily normalized to N(>Fx)=2.0x 10~^ at Fx=1.5xl0~^^ ergs~^cm~^. All three samples 
fall well below the Euclidean relationship at faint flux levels. The RGB catalog exhibits a higher 
density of objects than the HQS given identical flux limits. At a flux of ~10~^^ ergs~^cm~^, 
the HQS number density is ~1.7 times lower than the RGB sample likely a result of the very 
restrictive q;ox<1-1 (log(Sx/So)>1.3) HQS selection criterion (see Figure 3). 

The X-ray logN-logS distribution shows slight evidence for the "bump" at 
~3 X 10~^^ ergs~^cm~^ reported both by Maccacaro et al. (1989) and Nass et al. (1996), 
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Fig. 7. — The X-ray logN-logS relationship for the RGB Complete, EMSS and HQS samples 
of BLLacs. The latter two samples have been truncated to conform with the optical and radio 
selection criteria used in the creation of the RGB Complete sample. The solid line represents the 
no-evolution Euclidean NocS"^/^ powerlaw arbitrarily normalized. 

although the number statistics are poor. There also appears to be a rise in the surface density 
near Fx=8 x 10~^^ ergs~^cm""^. The number statistics are better here, but it is still difficult to 
assess the significance of this feature. The distribution could alternately be characterized by a 
powerlaw with a deficit of sources near 5 — 10 x 10~^^ ergs~^cm~^. 

Bade ct al. (1998b) have noted that the X-ray logN-logS distribution of HQS BLLacs is more 
monotonic when extreme HBLs (aox<0.91) are excluded. From this, Bade et al. (1998b) infer 
that intermediate BL Lacs constitute the "basic BL Lac population" and that X-ray dominated 
BL Lacs are those objects which are observed in a state of enhanced X-ray activity. In order to 
test this hypothesis, we examined the RGB X-ray logN-logS distribution, excluding first those 
RGB BLLacs with aox<0.91 as was done in Bade et al. (1998b). This criterion eliminates only 
two objects from the Complete sample, so the overall shape is not affected and the fcatTirc at 
Fx=8 X 10^^^ ergs~^cm~^ remains. Even if we raise the cutoff to values of Oqx to 1.0, 1.1 or 
1.2, the feature at Fx~8 x 10~^^ ergs~^cm~^ remains although the number statistics necessarily 
deteriorate. This contradicts the result of Bade et al. (1998b), although we cannot rule it out 
because of the unknown effects of the optical magnitude limit on the RGB Complete sample. 
Therefore we find the feature at Fx~8 x 10~^^ ergs^^cm"^ is significant, and does not seem to 
depend on properties of only an extreme subset of BL Lacs. We believe it therefore does not 
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Fig. 8. — The radio logN-logS relationship for the RGB Complete, EMSS, HQS and 1 Jy samples of 
BLLacs. The latter three samples have been truncated to conform with the selection criteria used 
in the creation of the RGB Complete sample. The solid line represents the no-evolution Euclidean 
N(>S)cxS^'^/^ powerlaw arbitrarily normalized. 

indicate extreme HBLs are objects in a state of enhanced X-ray activity, but suggest it is an 
indication of the type of breaks which are expected in logN-logS distributions of beamed objects 
(see below). 

5.6. Radio logN— logS Relationship 

Figure 8 shows a radio logN-logS distribution with the RGB, EMSS, HQS and 1 Jy samples 
included. As before, all samples have been truncated to include only those objects which adhere 
to the RGB selection criteria. This has the effect of reducing the 34 object 1 Jy complete sample 
to 19 objects. Note also that for all samples we plot core rather than total radio flux densities 
because the total emission contains kiloparsec-scale emission which is not as intimately connected 
with the relativistic flow in the center of the AGN. Many studies do not make this core vs. total 
flux distinction, but considering the wide range of core-to-lobe ratios exhibited by BL Lacs, we feel 
it important to do so here. 

Several points are evident from the figure. The incompleteness of the 1 Jy sample below radio 
core flux densities of ~1 Jy is apparent from the turnover in the number counts in the diagram. 



- 22 - 




10 ICQ 1000 

goore (mJy) 



Fig. 9. — The radio logN-logS relationship for the RGB Complete, EMSS and 1 Jy samples of 
BLLacs. The solid line represents the theoretical beaming model prediction of Urry, Padovani &; 
Stickel (1991; §5.6). 

However, above ~700mJy the HQS, RGB and (to within ~20%) the IJy samples agree well, 
indicating that the effects of the X-ray flux and optical magnitude limits do not do not severely 
affect the RGB Complete sample at these high radio flux density limits. In the middle part 

of the diagram (~70— 700 mJy), the RGB and HQS samples are roughly consistent with each 
other, unlike the results obtained for the X-ray logN-logS distribution. At the faintest radio flux 
densities, the surface density of HQS BL Lacs falls ~2 times below the RGB sample, indicating 
selection criteria based on identifying optically bright sources with flat values of Oqx produce an 
incompleteness that increases with decreasing radio flux density. 

While it may be tempting to infer the evolutionary behavior of BL Lacs from the curvature of 
the RGB logN-logS distributions, it is dangerous to do so without careful modeling of the effects 
introduced by the X-ray and optical magnitude flux limits. This is illustrated in Figure 9 which 
shows the radio logN-logS distribution for the RGB Complete sample in comparison with the full 
(not truncated) 1 Jy and complete EMSS samples. Here the effects of the radio, optical and X-ray 
limiting flux densities inherent in the RGB catalog are clearly evident, particularly at low radio 
flux densities where the counts flatten appreciably, falling a factor of ~4 below the EMSS sample, 
which itself represents only a fraction of the radio faint BL Lac population. 

Figure 9 also shows a theoretical beaming model for the logN-logS distribution (see and 
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Urry and Shafer 1984 and Urry et al. 1991 for details). This model is based on a simulation of a 

randomly oriented parent population of FR I radio galaxies which follow an assumed cosmological 
evolution model and whose jets have a powerlaw distribution of bulk jet Lorentz factors. The 
resulting BL Lac luminosity function is characterized by a double powerlaw, flattened by beaming 
at low radio powers with a break that is shifted to higher radio powers relative to breaks in 
the luminosity function of the parent population. When converted to surface number densities 
and fluxes, the logN-logS distribution of BLLacs exhibits a Euclidean slope (—3/2) down to 
100 mJy, flattens to N(>Sr)ocSjr"'^ between 10<Sr<100 mJy, and then becomes very flat (ocS~'^'^) 
for Sr<10mJy. Beaming therefore tends to flatten source counts to an increasing degree at lower 
flux densities. 

Unfortunately, because neither the EMSS nor RGB samples are truly complete (containing 
all BLLacs) at radio flux densities below ~lJy, it is not possible to directly compare the specific 
predictions of this model with any of the samples in Figure 9 other than the 1 Jy which provides 
the normalization. However, we can make some general observations. The RGB sample follows 
a Euclidean slope above ~200 mJy. At intermediate radio flux densities, it flattens more than 
the model but, as stated above, this is at least partly a result of the incompleteness caused by 
the RGB Complete sample's flux limits. At ~40mJy, the surface density of the RGB and EMSS 
samples coincide, but below ~20 mJy, the Complete EMSS (HBL) sample quickly overwhelms the 
RGB sample. Nevertheless, the EMSS sample falls well below the prediction of Urry et al. (1991). 
It is this deficit of low flux/luminosity BLLacs which led Morris et al. (1991) to postulate that 
HBLs may evolve negatively, a claim reinforced by Rector et al. (1999a) in their new analysis of 
the V/V 

max analysis of an updated Complete EMSS sample. However, the EMSS survey consists 
of only the most extreme X-ray— dominated objects. As discussed below, we believe the dichotomy 
of HBLs and LBLs is a result of selection effects and conclusions drawn about the evolution of 
only one extreme end of the distribution should therefore be treated cautiously. 



6. Discussion 

6.1. Comparison with Previous BLLac Surveys 

The principal result of this study is that the RGB BL Lacs exhibit a continuous range in SEDs 
rather than segregating into two distinct classes: HBLs vs. LBLs (or XBLs vs. RBLs). Clearly the 
RGB sample smoothly spans a vast range of SEDs: 5 orders of magnitude in Sx/Sr, 0.5<aox<l-7 
and 0.2<aro<0.8. Most RGB BLLacs have intermediate properties (e.g., Sx/Sr ratios) between 
the HBL and LBL extremes that dominated earlier samples. Similar preliminary results have also 
been reported for the deep DXRBS and REX X-ray surveys (Perlman et al. 1998; Caccianiga et al. 
1999). The question remains whether these results are an accurate representation of the BLLac 
population as a whole, or whether the true BL Lac distribution is bimodal and the flux limits 
of the RGB survey unfortuitously tuned to make make it particularly sensitive to objects with 
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Fig. 10. — Comparison of the log Sx/Sr ratio for the full RGB sample (upper solid line) and those 
objects with a. radio flux densities greater than 500 mJy (hatched histogram) or b. X-ray fluxes 
greater than 3.5 x 10~^^ ergs~^cm~^. Note the similarity between the hatched distributions and 
those in Figure 4. 

intermediate SEDs. Occam's Razor makes it tempting to hypothesize that the intermediate nature 
of the RGB sample is an accurate reflection of an underlying unimodal distribution and that the 
previously observed bimodal distribution is a result of high flux limits of previous surveys, but 
this hypothesis bears closer examination. 

Figure 10 shows what the RGB's distribution of Sx/Sr ratios would have been had either the 
radio (Fig. 10a) or X-ray (Fig. 10b) flux limits been higher. Not surprisingly, increasing the radio 
flux limit preferentially selects increasingly more radio-dominant BL Lacs. At a limiting radio flux 
density of 500 mJy a distribution roughly consistent with that of the 1 Jy sample (Fig. 4) results. 
(Too few RGB objects arc brighter than 1 Jy to make a meaningful comparison at this higher radio 
flux density limit.) Direct comparisons with X-ray— selected samples are complicated since none of 
the comparison surveys were conducted over the same energy band as the RGB (0.1— 2.4keV with 
the ROSAT PSPC) and X-ray flux limits are sensitive to the assumed photon index and Galactic 
column density. In addition, the Slew Survey has a factor of several range of limiting sensitivities 
while the HEAO-1 sample is a result of a flux density (not flux) limited survey complicating a 
precise comparison with the RGB sample. 

Despite these ambiguities, increasing the limiting X-ray flux of the RGB survey selects 
the more X-ray— dominated RGB BLLacs and at a flux limit of 3.5 x 10^^^ ergs^^cm^^, the 
distribution is roughly consistent with both the HEAO-1 and Slew Survey distributions (Figs. 4a 
and 4b). However, the effective limit of the EMSS in the ROSAT band is 2.2 - 3.5 x lO^^ ^jy (T. 
Rector, private communication) which is only marginally brighter than the RGB X-ray flux limit. 
Therefore if our assumption of a unimodal Sx/Sr ratio for BLLacs is correct, then we predict 
(based on only a small difference in the RGB and EMSS X-ray survey flux limits), that the EMSS 
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sample should closely resemble that of the RGB; this clearly is not the case. However, Rector 

et al. (1999a) have made a careful re-examination of the EMSS identifications and have found 
several new BL Lacs that were misidentified because of their strong 4000 A break contrasts. The 
addition of these sources does bring the EMSS and RGB Sx/Sr distributions into somewhat closer 
agreement. Also, a new analysis by Rector et al. (1999b) which extends the complete EMSS 
sample down to 2.0 x 10~^^ ergs~^cm~^ and declinations >40°, also preferentially adds objects 
with intermediate Sx/Sr ratios, but not sufficiently to make the RGB and EMSS distributions 
consistent. Nevertheless, agreement between the 1 Jy, HEAO-1 and Slew samples with that of the 
RGB sample were either its radio or, respectively, X-ray flux limits raised, is compelling. 

Another demonstration that the RGB survey likely accurately reflects an underlying unimodal 
distribution arises from the following Monte Carlo simulation. We create hypothetical populations 
of BL Lacs where both radio and X-ray log A^— log S distributions obey simple powerlaws, and radio 
and X-ray fluxes are independent variables. The first assumption is clearly an oversimplification 
since there is ample evidence for breaks in both the radio and X-ray log AT— log S relationships 
(including this paper) but including these breaks would introduce many new parameters to the 
model (location of the break(s) and the slopes above an below them) without adding substantially 
to the general sense of the results. Figure 11 shows a simulation of 30,000 sources whose radio and 
X-ray fluxes are constrained to be >lmJy and >5 x 10~^^ ergs~^cm~^, respectively. These limits 
were chosen because they are fainter than either the radio or X-ray flux of any BL Lac in any of 
the relevant comparison samples (EMSS, RGB or IJy). We then applied the different flux limits 
of the comparison surveys and randomly chose objects in order to create samples with the same 
number of objects. The particular simulation shown assumed the log A— log S powcrlaw slopes 
were —1.2 and —0.8 in the radio and X-ray, respectively. (See Urry, Padovani &; Stickel 1991 and 
Maccacaro et al. 1984 for the relevance of these specific indices.) 

Clearly Figure 11 suggests the underlying distribution is unimodal; this is, in fact, a universal 
among all the simulations, although the location of the peak in the distributions varies with the 
radio and X-ray logN— logS slopes. The distribution shown in Figure 11 peaks at logSx/Sr=— 5.61. 
The distributions resulting from flux limits at IJy (5 x 10"^^ ergs~^cm~^) are acceptably 
consistent with the distributions shown in Figure 4 for the 1 Jy (EMSS) samples while the sample 
created by enforcing a radio flux limit of 20mJy and an X-ray flux limit of 3 x 10~^^ ergs^^cm^^ 
is an acceptable match to the RGB sample's distribution, given the simplicity of our assumptions. 

Both the above analyses strongly point towards the simple conclusion that the true distribution 
of BLLac SEDs, and Sx/Sr ratios in particular, is unimodal and accurately represented by the 
RGB survey. Our results also show that the previously observed bimodality can be explained as 
a result of naturally occurring observational selection effects present in older surveys. A similar 
conclusion is reached by Caccianiga et al. (1999) in their preliminary examination of the REX 
survey BLLacs and in the theoretical modeling of Fossati et al. (1997). 
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Fig. 11. — Results of a Monte Carlo simulation of the logSx/Sr for BLLacs. The model shown 
assumes the log iV— log S relationship follows a S"-*^-^ powerlaw in the radio and a S^'^-^ powerlaw in 
the X-ray. a. The underlying distribution which is unimodal and peaks at log Sx/Si=— 5.6 b. The 
distributions which would result from this sample if surveys with flux limits of 1 Jy (right-slanted 
hatched histogram), 5 x 10"^"^ ergs~^cm~^ (left-slanted hatched histogram) or dual flux limits of 
20 mJy and 3 x 10~^^ ergs~^cm~^ were performed. Flux limits have been chosen based on the limits 
inherent in the 1 Jy, EMSS and RGB surveys, respectively and the total number of objects in these 
simulated samples have been matched to the number of objects in the survey they are intended to 
represent. Although the agreement is not perfect, the similarity with Figure 4 is suggestive that 
the logSx/Sr distribution for BLLacs is unimodal and the bimodal distribution reflected in the 
EMSS and 1 Jy samples is a result of the flux limits of these different surveys. 



6.2. A Unified BL Lac Population 

Differences between the various subclasses of BLLacs have been attributed to different jet 
orientations to the line-of-sight (e.g., Kollgaard 1994; Urry and Padovani 1995). However, further 
analysis showed that the range of SEDs observed in LBLs and HBLs cannot be reproduced by 
simple changes in jet orientation (e.g., Sambruna et al. 1996; Fossati et al. 1997) and our results 
confirm this (§5). Padovani and Giommi (1995a) propose an alternative to this orientation-based 
model, namely, that LBLs and HBLs essentially share the same range in orientation, but have 
intrinsically different SEDs. They postulate that the frequency at which the synchrotron break 
occurs, t^rcak) differs intrinsically in LBLs and HBLs. This then alters the balance of synchrotron 
and inverse Compton emission present in soft X-rays, moving the objects in the Oro— ctox diagram. 
This model predicts that BL Lacs whose spectra break at high frequencies should be HBL-like and 
lie along the horizontal path labeled "HBLs" in Figure 3. Objects with i^reak at low frequencies 
should lie along the diagonal swath labeled "LBLs" in the aro— c^ox plane. Assuming (ctox)— 0.7 and 
(ckro)— 0.37, as derived from the EMSS sample, and a spectral index abrcak~l-9, valid for z^>i4)rcak) 
Padovani and Giommi (1995a) were able to reproduce the bimodal locations of the EMSS XBLs 
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and IJy RBLs in the color-color diagram. If this model is correct, then all new BLLacs should 

continue to reside along constrained loci in the ajo— ckox plane, although this SED-bascd model 
need not lead to a bimodal separations of classes. Our results show that this may be the case with 
the RGB objects following the traditional paths in the a^o vs. aox plane, although the scatter is 
significant and the triple-flux limits restrict the range of SEDs to which we are sensitive. 

Two additional new unification models attempt to explain the physical origin of the wide 
range of SEDs in BL Lacs. The first asserts that the shape of the SED is linked to the bolometric 
luminosity (Fossati et al. 1997). In this model, HBLs and LBLs are manifestations of the same 
phenomenon, but the broadband SED varies in a predictable way with bolometric luminosity. 
This model can reproduce the disparate distribution of the 1 Jy and Slew Survey samples in the 
CKro~ctox plane, and predicts an intrinsically smooth transition between the two extremes. The 
other model predicts that the origin of the SED differences lies in differences in the electron kinetic 
luminosity of the jet which is related to jet size (Ljet oc r^; Georganopoulos and Marscher 1998). 
This model also predicts a smooth transition between the previously disparate BL Lac subclasses 
which qualitatively agrees with the distribution of RGB objects. Quantitative comparison of these 
models with our data is, however, difficult both because the RGB sample is triply flux limited 
and because the RGB Complete sample's limiting magnitude of 18.0 curtails the range of possible 
SEDs to which the sample is sensitive. We are however currently attempting to model the selection 
effects, but extending the Complete sample to a fainter magnitude limit would clearly alleviate 
some restrictions. 

The wide range of spectral indices exhibited by the RGB BL Lacs has important implications 
for followup surveys for RASS BL Lacs. Figure 3 shows not only the distribution of BL Lacs in 
the CKro— CKox diagram, but also the color-color classification boundaries established in the EMSS 
by Stocke et al. (1991). These regions have been used to successfully identify new BLLac objects. 
Nearly all Einstein Slew Survey objects in the region defined by aio=(0.3,0.6) and aox=(0.55,1.2) 
were spectroscopically confirmed as BLLacs (Schachter et al. 1993b; Perlman et al. 1996a). A 
criterion of q;ox<1-1 was used to select candidates for the Hamburg Quasar Survey BLLac sample 
(Nass et al. 1996). However, the RGB sample, selected without any spectral index criteria, shows 
that targeting only those objects in the RASS with particular color-color indices will miss a large 
fraction of the BLLacs. If, for example, the EMSS BLLac class boundaries had been used to 
select candidates, 67% of the RGB BLLacs would have been excluded (70% of the sources in our 
Complete sample). Using the criterion in Nass et al. (1996) would have excluded 70% of our 
objects (and 76% of objects in the Complete sample). Therefore while the likelihood of finding 
BLLacs does increase dramatically with decreasing aox, candidates cannot be chosen based on 
their SEDs without producing highly biased samples leading to incorrect or uncertain conclusions. 

Finally, because we lack detailed knowledge of the effects of the triple flux limits on the 
characteristics of the entire RGB BL Lac sample, we are not at present able to definitively evaluate 
the predictions of any of the FR I/BLLac or HBL/LBL unification models. Some insights might 
be obtained with Monte Carlo simulations which exclude sources which fall below our RGB fiux 
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limits, but this, of course, requires a more elaborate simulation of the intrinsic properties of the 
BLLac population than the simple simulation shown in Figure 11. Two additional approaches 
can be pursued. First, redshifts for the RGB Complete Sample can be obtained with currently 
available optical telescopes. This would transform the source counts into correctly normalized 
luminosity functions, which could then be compared with luminosity functions of FRI radio 
galaxies. But here again, the selection effects must be precisely taken into account. Second, the 
RGB sample can be extended to reduce the effects of the flux limits. Most importantly, spectra 
should be obtained for RGB objects fainter than ~18.5 in order to discriminate BLLacs from 
other X-ray/radio sources. While improvements in the X-ray flux limit the RGB catalog are not 
likely in the near future, the NVSS (Condon et al. 1998) and FIRST (Becker et al. 1995) radio 
surveys permit a factor of 2±20 reduction in the radio flux limit. BL Lac samples based on these 
radio surveys are currently being pursued (Caccianiga et al. 1999; Laurent- Muehleisen et al. 1999). 
Here again, substantial spectroscopy of optically faint objects will be needed. Significant progress 
in refining BLLac unification scenarios thus depends on spectroscopy with new 8-meter class 
telescopes, which fortunately are proliferating today. 
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09 


1 fi 

10 


c;i Q9 
01. yz 


KO 
OZ 


^8 
oo 


98 

zo. 


.4 


46 


5.35 


17A 


1 

1 , 


1 8 

. lo 


0.37 


1 


.14 


0.190 


N96,L97,B98 i 


0928+747 


8C 0923+750 


09 


98 


09 07 


7A 

I 4 


A7 

4 / 


1 


.1 


5 


1.98 


21.0 


n 

u. 


78 
. ) o 


0.43 


0, 


.72 


0.638 


B95,P96,L97 


0929+502"^ 


RXJ09292+5013 


09 




1 5 AA 


^n 




ou. 


.0 


895 


1.19 


17.2 


1 


^n 

.ou 


0.60 


1 


.43 




N96,B97,L97 


0930+350*^ 


7C 0927+3516 


09 


30 


55.28 


35 


03 


37, 


.6 


337 


0.71 


20.4 






0.77 


1 


.03 




H97,L97 


0952+656 


RXJ09524+6538 


09 


52 


32.19 


65 


38 


01, 


.1 


27 


1.41 


18.0 


1 


.23 


0.37 


1 


.27 




L98,B97,L97 


1000+225A 


4C+22.25 


10 


00 


21.36 


22 


33 


07, 


.4 


<5 


0.96 


17.8 






<0.20 


1 


.35 


0.419 


B95,L97 


1012+424'^ 


B3 1009+427 


10 


12 


44.30 


42 


29 


57, 


.2 


29 


7.68 


18.1 


0, 


.84 


0.39 


0, 


.92 




L98,B97,L97 


1015+494 


IH 1011+496 


10 


15 


04.17 


49 


26 


00, 


.6 


242 


19.5 


16.5 


1 


.49 


0.44 


1 


.14 


0.2? 


B95,L96,L97 


1047+546*^ 


lES 1044+549 


10 


47 


45.82 


54 


37 


41, 


.3 


6 


0.45 


21.2 


0, 


.63 


0.49 


0, 


.89 




PG95,B97,L97 


1053+494*^ 


MS 1050.7+4946 


10 


53 


44.14 


49 


29 


56, 


.1 


19 


1.20 


17.5^ 


1, 


.62 


0.30 


1, 


.46 


0.140 


R99,B97,L97 


1104+382 


MKN 421 


11 


04 


27.31 


38 


12 


31, 


.8 


600 


405. 


13.3 


1, 


.37 


0.28 


1, 


.11 


0.031 


B95,L96,VLA 


1120+422 


lES 1118+424 


11 


20 


48.07 


42 


12 


12, 


.5 


19 


6.65 


17.1 


1 


.30 


0.27 


1 


.18 


0.124? 


PG95,L98,B97,L97 


1149+246^* 


EXO 1149.9+2455 


11 


49 


30.3 


24 


39 


27 




15 


9.12 


18.9 






0.39 


0, 


.85 




PG95,L97 


1209+413 


B3 1206+416 


12 


09 


22.79 


41 


19 


41, 


.4 


393 


0.87 


17.2 


1 


.10 


0.53 


1 


.43 




H97,L97 



TABLE 2— Continued 



RGB J 
Name 


Alternate 
Name 




RA 

(J2000) 


Dec 




Qcore 

(mJy) 


Fx 



(mag) 




Q!ro 


a 


ox 


z 




Ref. 


1215+075*^ 


lES 1212+078 


12 


10 


1 n Q 


07 


'^9 
oz 


04 
U4 




84 


4.39 


I 7 

I I , 


8 

.o 


n 

u. 


88 
.oo 


0.45 


1. 


.06 


0. 


,130 


PG95,B97,L97 


1217+301 


lES 1215+303 


12 


1 7 


f^9 ns 

OZ.Uo 


oU 


07 
U ) 


00 

UU. 


,6 


353 


38.4 


1 

10, 


.4 


1 

1 , 


88 
.oo 


0.38 


1 


.29 


0. 


,237 


B95,L96,L97 


1220+345 


RX J12201+3431 


12 


on 


ns ^9 


^4 


^1 
ol 


91 
Z 1 . 


,6 


258 


0.42 


1 8 

lo, 


R 
.D 


n 

u. 


Q8 
.yo 


0.61 


1 


.32 






B95,L96,L97 


1221+301^^ 


lES 1218+304 


12 


91 


91 Q4 


OU 


1 n 

lU 


'^7 
O 1 . 


1 


60 


21.6 


1 fi 

ID, 


7 


1 

1 , 


99 
.zz 


0.34 


1 


.04 


0. 


,182^ 


B95,L96,L97 


1221+282 


W COMAE 


12 


91 


^1 fiQ 


98 


1 ^ 
lo 


c;8 
oo. 


,5 


940 


2.25 


1 

10, 


oa 
.o 


1 

1 , 


94 

.Z4 


0.50 


1 


.51 


0. 


,102 


B95,L96,K92 


1224+246^* 


MS 1221.8+2452 


12 


94 


94 9 


94 


oD 


94 

Z4 




25 


2.38 


1 

ID, 


Q 

.y 


1 

1 , 


47 
.4 ( 


0.28 


1 


.40 


0. 


,218 


B95,P96,L97 


1230+253 


RXJ12302+2517 


12 


OU 


14 07 
14. U 1 


9c; 

ZO 


1 8 
10 


07 
U ) . 


,0 


351 


3.35 


I 7 

I I , 


.4 


1 

1 , 


8Q 
.oy 


0.54 


1 


.40 




c 
I 


N96,B94,L97 


1310+323 


AUCVn 


13 


1 n 

iU 


98 fifi 
Zo.DD 


'^9 


90 
ZU 


4*? 
40. 


,8 


1970 


0.69 


1 Q 

ly , 


Q 

.y 


n 

u. 


.yo 


0.83 


1 


.06 


0. 


,996 


B95,U96,G93 


1341+399'^'^ 


RXJ13410+3959 


13 


41 


uo. lo 


«3y 


oy 


40. 


,5 


34 


7.92 


1 8 

lo. 


9 
.z 






0.41 


0, 


.97 


0. 


,163 


N96,L97 


1417+257*^ 


EXO 1415.6+2557 


14 


1 7 
1 ) 


00. D 


9^; 
zo 


4*? 


9f^ 
ZO 




40 


23.7 


1 7 

1 ) , 


.0 


1 

1 , 


1 Q 

. ly 


0.36 


0, 


.90 


0. 


,237 


B95,L98,L96,L97 


1422+580 


lES 1421+582 


14 


99 
zz 


■^s on 


c;8 

Oo 


01 
Ul 


00. 


,5 


6 


22.0 


1 8 

lo. 


7 


1 

1 , 


04 
.U4 


0.27 


0, 


.71 


0. 


,638 


B95,B94,L97,B98 


1427+238 


PG 1424+240 


14 


97 


nn "io 
uu.oy 


9*? 
Zo 


48 
4o 


00 

uu. 


,0 


250 


6.66 


1 fi 

ID, 


/I 
.4 


1 

1 , 


48 
.4o 


0.43 


1 


.32 






B95,K96,L97 


1439+395 


PG 1437+398 


14 


oy 


1 7 4Q 
1 ) .4:y 


oy 


'?9 
oz 


49 

4Z. 


,9 


38 


15.3 


1 8 

lo. 


1 

. 1 


1 

1 , 




0.41 


0, 


.97 






B95,L98,F93,L97 


1442+120*^ 


lES 1440+122 


14 


49 


48 ^ 


1 9 


no 


40 

4U 




45 


11.0 


1 7 

1 1 , 


4 


1 


90 


0.37 


1 


.04 


0. 


,162 


B95,L98,L97 


1448+361 


RXJ14479+3608 


14 


48 


00.58 


36 


08 


31. 


,7 


29 


6.83 


17, 


.2 


1 


.60 


0.32 


1 


.23 






L98,B97,L97 


1454+514*^ 


RX J14544+5124 


14 


54 


27.12 


51 


24 


33. 


,6 


80 


0.95 


18, 


.4 






0.49 


1 


.27 






L98,L97 


1509+559 


RX J15098+5555 


15 


09 


47.96 


55 


56 


17. 


,4 


28 


1.08 


18, 


.4 


1 


.99 


0.41 


1 


.46 






K96,B97,L97 


1534+372^^ 


RXJ15347+3716 


15 


34 


47.20 


37 


15 


54. 


,8 


20 


0.49 


18, 


.3 


1 


.84 


0.37 


1, 


.55 


0. 


,143 


L98,B97,L97 


1536+016^^ 


MS 1534.2+0148 


15 


36 


46.8 


01 


37 


59 




35 


8.17 


19 


.9 


0, 


.89 


0.53 


0, 


.67 


0. 


,312 


B95,P96,L97 


1540+147 


4C+14.60 


15 


40 


49.49 


14 


47 


45. 


,9 


716 


1.90 


16 


.4 


0, 


.66 


0.50 


1 


.36 


0. 


,605 


B95,U96,L97 


1544+049^* 


RXJ15442+0458 


15 


44 


18.7 


04 


58 


22 




39 


1.70 


19, 


.4 


0, 


.74 


0.51 


0, 


.95 






K96,B97,L97 


1555+111 


lES 1553+113 


15 


55 


43.05 


11 


11 


24. 


,4 


398 


38.7 


15 


.5'^ 


1, 


.45 


0.39 


1, 


.16 


0. 


,360 


B95,K96,L97 


1602+308 


RXJ16022+3050 


16 


02 


18.09 


30 


51 


09. 


,3 


20 


2.11 


18, 


.5'' 


1 


.56 


0.38 


1 


.22 






L98,K96,L97 


1621+377 


4C+37.46 


16 


21 


11.29 


37 


46 


04. 


,9 


155 


1.23 


21, 


.4 






0.78 


0, 


.81 






K96,L97 


1624+374 


B3 1622+375 


16 


24 


43.35 


37 


26 


42. 


,4 


14 


0.25 


19 


.4 






0.42 


1 


.33 


0. 


,200 


K96,L97 



TABLE 2— Continued 



RGB J 


Alternate 




RA 




Dec 


ccore 
'-'5 


Fx 







Q!ro 


"ox 


z 


Ref. 


Name 


Name 






(J2000) 






(mJy) 




(mag) 












1626+352 


RX J16264+3513 


16 


26 






13 


41.5 


14 


1.16 


18.7 


0.89 


0.35 


1.14 


0.500 


K96,K96,L97 


1644+457 


B3 1642+458 


16 


44 


1Q QS 


4"^ 


46 


44.5 


64 


2.60 


18.3 


1.14 


0.46 


1.12 


0.220 


N96,K96,L97,B98 


1719+177 


PKS J1719+1745 


17 


19 




1 7 


45 


06.4 


600 


0.87 


19.2 




0.73 


1.17 




B95,VLA 


1725+118 


IH 1720+117 


17 


25 


04. ^8 

Ut:.00 


1 1 


52 


15.5 


88 


36.2 


16.3" 




0.35 


1.01 


0.018 


B95,L97 


1743+195 


lES 1741+196 


17 


43 


^7 84 
O 1 .o^ 


1 Q 


35 


08.9 


157 


9.69 


17.3" 


0.98 


0.47 


1.03 


0.083 


PG95,B97,L97 


1747+469 


B3 1746+470 


17 


47 




4fi 


58 


50.9 


548 


0.40 


21.3 


2.08 


0.82 


1.10 


>1.484 


V96,L97 


1750+470^ 


B3 1748+470 


17 


50 


0^ fin 


47 


00 


43.7 


10 


4.01 


19.0^^ 


0.94 


0.37 


0.91 


0.160 


L98,B97,L97 


1808+468 


RX J18079+4648 


18 


08 




4fi 


49 


41.0 


41 


0.93 


18.4 


1.70 


0.44 


1.40 




L98,L97 


1811+442 


RX J18118+4416 


18 


11 


<^'i 47 


44 


16 


28.5 


6 


0.64 


18.9" 


1.18 


0.31 


1.26 




L98,B97,L97 


1813+317 


EXO 1811.7+3143 


18 


13 


35.21 


31 


44 


17.7 


74 


0.99 


17.9" 




0.45 


1.33 


0.117 


PG95,L97 


1824+568 


4C+56.27 


18 


24 


07.07 


56 


51 


01.5 


1120 


2.76 


18.2 


0.15 


0.67 


1.02 


0.664 


B95,U96,K92 , 


1829+540 


RX J18293+5403 


18 


29 


24.29 


54 


02 


59.8 


18 


2.35 


17.6 


2.31 


0.31 


1.55 




L98,B97,L97 w 


1838+480 


RX J18387+4802 


18 


38 


49.17 


48 


02 


34.4 


23 


2.78 


17.6 


2.10 


0.33 


1.46 




L98,L97 1 


1841+591 


RXJ18412+5906 


18 


41 


20.31 


59 


06 


08.2 


6 


0.74 


19.7^ 


2.54 


0.35 


1.44 


0.530 


L98,B97,L97 


1848+427 


RXJ18487+4245 


18 


48 


47.14 


42 


45 


39.4 


8 


12.2 


18.9 


1.20 


0.34 


0.81 




L98,L97 


2039+523"^ 


lES 2037+521 


20 


39 


23.50 


52 


19 


49.9 


19 


4.06 


19.5" 


1.80 


0.46 


1.05 




PG95,B97,L97 


2145+073 


MS 2143.4+0704 


21 


45 


52.30 


07 


19 


27.2 


37 


8.00 


19.0*^ 


1.91 


0.47 


1.04 


0.237 


B95,L98,P96,L97 


2202+422 


BL LACERTAE 


22 


02 


43.29 


42 


16 


40.0 


2308 


5.68 


15.1" 


0.94 


0.53 


1.41 


0.0686 


B95,U96,L97 


2250+384^ 


B3 2247+381 


22 


50 


05.77 


38 


24 


37.3 


60 


7.36 


16.0"'*^ 


1.51 


0.29 


1.37 


0.119 


L98,B97,L97 


2322+346^^ 


RXJ23226+3436 


23 


22 


44.01 


34 


36 


14.0 


30 


2.83 


18.2^= 


0.80 


0.40 


1.05 


0.098 


L98,B97,L97 


The Complete Sample 


0710+591 


IH 0658+595 


07 10 30.07 


59 08 20.5 


34 


25.2 


18.0^= 


1.10 


0.39 


0.81 


0.125 


PG95,L98,B97,L97 


0721+713 


IH 0717+714 


07 21 53.45 


71 20 36.4 


600 


3.68 


15.4 


2.02 


0.43 


1.70 




B95,B97,VLA 


0806+595 


RXJ08063+5931 


08 


06 25.94 


59 31 06.9 


29 


4.26 


17.2 


0.50 


0.32 


1.09 




L98,B97,L97 


0958+655 


S4 0954+65 


09 58 47.24 


65 33 54.8 


480 


1.10 


16.7 


0.24 


0.50 


1.35 


0.367 


B95,U96,G93 



TABLE 2— Continued 



RGB J 

Name 


Alternate 
Name 




RA 

(J2000) 


Dec 




Qcore 
^5 

(mJy) 


Fx 

(10-12^) 



(mag) 




"ro 


a 


OX 


z 




Ref. 


1031+508 


lEb 1028+511 


10 


oi 


1 s c;9 


OU 


Oo 


00. 


,8 


23 


r A o 

59.2 


16.3 


1 

1 . 


^7 
.0 ) 


O O I 

0.21 


0, 


A/* 

.96 


0. 


O/^ 1 

.361 


TiAr T AO "DAT T AT "nAO 

B95,L98,B97,L97,B98 


1037+571 


KX J10377+5711 


10 


O ) 


4zl "il 


t 


1 1 
i 1 


00. 


,8 


89 


o r o 

2.58 


t I"? A 


1 

1 , 


8n 

.oU 


O /I o 

0.43 


1 


A 1 

.41 






ATA/^ T AO "OAT T AT 

iN9D,L98,B97,L97 


1136+701 


MKN 180 


11 


OD 


9fi 41 

ZU.4:1 


7n 


no 
uy 


97 
Z 1 . 


,3 


139 


58.3 


15.5 


1 

1 . 


1^1 

.01 


O OO 

0.32 


1 


I O 

.13 


0, 


.046 


TiAr" T A/^ T AT 

B95,L96,L97 


1136+676 


RX J11365+6737 


11 


ou 


ou.uo 


fi7 


■^7 


04 
u^. 


,4 


A r\ 

40 


24.8 


17.4 


1 

1 , 


9Q 

.zy 


0.36 


0, 


.93 


0, 


.136 


T AO T~i f\ A T AT T>AO 

L98,B94,L97,B98 


1151+589 


8C 1148+592 


11 


01 


9zL fiS 


oo 


oy 


1 7 
1 ) . 


,6 


95 


1.66 


17.8 


n 

u, 


fi9 
.DZ 


0.46 


1, 


.16 






TVTA/^ nAT T AT 

N96,B97,L97 


1231+642 


Mb 1229.2+6430 


12 


1 


1 .^u 


fi4 


14 
1^ 


1 8 

lo. 


,4 


34 


4.33 


17.0 


n 

u, 


QQ 

.yy 


0.31 


1 


.19 


0, 


.170 


"I~> A 1^ TIA/^ T AT 

B95,P96,L97 


1248+583 


PG 1246+586 


12 




lo. 1 o 


c;8 

Oo 


9n 
zu 


98 
Zo. 


,7 


204 


o or 

3.25 


■if A 

15.4 


1 

1 . 


49 

.4:Z 


O O /I 

0.34 


1 


f f 

.55 






TiAr T AT nAO 

B95,L97,B98 


1415+485 


RX J14155+4830 


1 A 

14 




89 


48 


ou 


ou. 


,5 


58 


r\ A A 

0.44 


17.6 


1 

1 , 


4n 

.4U 


0.40 


1 


.55 






T AO noT T AT 

L98,B97,L97 


1419+543 


b4 1418+54 


14 


1 Q 




0^ 


9*? 
zo 


14 


,8 


■1 A r\r\ 

1400 


-t r»o 

1.22 


■i f r\ 

15.0 


1 

1 . 


1 9 

. iz 


O /IT 

0.47 


1 


A 

.69 


0, 


.151 


B95,U96,VLA 


1427+541 


RX J 14274+5409 


1 A 

14 


97 


ou.oo 




uy 


9*? 
zo. 


,5 


O A 

24 


0.80 


"IT AC 

17 A 


n 

u, 


8fi 
.ou 


0.32 


1, 


O T 

.37 


0, 


.105 


T AO "r>A'V T A'V 

L98,B97,L97 , 


1428+426 


IH 1430+423 


14 


zo 


^9 fi9 


49 


4n 
^u 


91 

Zl. 


1 


O 1 

21 


61.4 


■i O C\C 

16.9 


1 

1 . 


.uo 


C\ 0^7 

0.27 


0, 


.81 


0, 


.130 


TiAr" T AO T A/^ T AT 

B95,L98,L96,L97 ^ 


1516+293 


T> "V TIP'!/'/' 1 

RX J15166+2917 


15 


1 fi 
10 


^i.oy 


9Q 


1 8 

lo 


no 
uy . 


,5 


34 


1.94 


18.2*^ 


1 

1 . 


'^9 
.oZ 


0.41 


1 


.22 


0, 


.130 


T AO noT T AT 

L98,B97,L97 i 


1517+6540 


IH 1515+ddO 


15 


1 7 
1 1 


47 fi8 




9^; 

zo 


9"? 
zo. 


,3 


19 


18.8 


17.4 


1 

1 , 


1 fi 
. lU 


0.26 


0, 


A /I 

.94 


>0.70 


T~> f\ f T AO nAT T AT "D A A 

PG95,L98,B97,L97,B99 


1532+302" 


RX J15319+3016 


15 


^9 


09 94 


ou 


Ifi 


9Q 
^y . 


,0 


A T 

47 


f o^y 

5.87 


1 r f c 

15.5 


1 


00 
. uu 


f\ OO 

0.23 


-< 

1, 


O/^ 

.36 


0, 


.064 


T A O TAT 

L98,L97 


1533+342 


RX J15333+3416 


15 


33 


24.25 


34 


16 


40. 


,2 


33 


5.49 


17.9 


1 


.60 


0.38 


1 


.17 






T AO TIAT T AT 

L98,B97,L97 


1542+614 


RX J 15428+6129 


15 


42 


56.94 


61 


29 


55. 


3 


102 


1.46 


1 T O 

17.2 


1 


.50 


O yl O 

0.42 


1 


A A 

.44 






T AO nAT T AT 

L98,B97,L97 


1554+/01 


Mb lo5/.l+/020 


15 


54 


24.13 


20 


11 


25. 


,4 


OO 

zO 


5.53 


17.0 


0, 


.79 


O OO 

0.30 


1, 


O CT 

.05 


0, 


.222 


noCT T)A£? T AT 

B95,P96,L97 


1652+403 


RXJ16527+4023 


16 


52 


49.93 


40 


23 


09. 


,9 


11 


0.49 


17.3 


1 


.80 


0.24 


1 


.68 






L98,L97 


1653+397 


MKN 501 


16 


53 


52.22 


39 


45 


36. 


,6 


450 


72.8 


14.2 


1 


.35 


0.32 


1 


.24 


0, 


.033 


B95,U96,G93 


1704+716 


RXJ17047+7138 


17 


04 


46.98 


71 


38 


17. 


,6 


17 


4.28 


16.5 


2, 


.10 


0.22 


1 


.55 






N96,L98,L97 


1728+502 


IH 1727+502 


17 


28 


18.62 


50 


13 


10. 


,5 


168 


31.3 


15.2 


1, 


.39 


0.32 


1, 


.24 


0, 


.055 


B95,L96,L97 


1742+597 


RXJ17424+5945 


17 


42 


32.00 


59 


45 


06. 


,8 


77 


0.64 


VIA 


1 


.96 


0.41 


1, 


.67 






L98,B97,L97 


1745+398 


B3 1743+398B 


17 


45 


37.76 


39 


51 


30. 


,8 


118 


1.59 


18.3*^ 






0.51 


1 


.20 


0, 


.267 


L98,L97 


1748+700 


S5 1749+70 


17 


48 


32.84 


70 


05 


50. 


,8 


610 


2.14 


16.9 


1 


.77 


0.51 


1. 


,45 


0. 


,770 


B95,U96,K92 


1749+433 


B3 1747+433 


17 


49 


00.36 


43 


21 


51. 


3 


281 


0.47 


17.9 


1 


.32 


0.56 


1. 


,47 






V96,L97 



TABLE 2— Continued 



RGB J 


Alternate 


RA Dec 


Qcore 
'-'5 


Fx 











Ref. 


Name 


Name 


(J2000) 


(mJy) 


(10-12^) 

\ s cm-' ' 


(mag) 










1756+553 


RXJ17561+5522 


17 56 15.89 55 22 18.1 


10 


13.6 


17.9 


1.64 


0.28 


1.04 ••• 


L98,B97,L97 


1806+698 


3C371 


18 06 50.68 69 49 28.1 


950 


2.51 


14.9° 


1.28 


0.44 


1.64 0.0512 


B95,U96,G93 


2243+203 


RXJ22438+2020 


22 43 54.73 20 21 03.9 


86 


0.97 


16.0 


1.68 


0.31 


1.72 ••• 


L98,B97,L97 


2319+161 


RXJ23196+1611 


23 19 43.44 16 11 50.1 


17 


3.23 


17.9" 


1.20 


0.32 


1.15 ••• 


L98,B97,L97 



^Optical magnitude converted from another band to equivalent O mag (aopt=1.0); See notes below 
^ID Uncertain; see notes below 

•^Optical magnitude corrected to include AGN component only 
'^Prom low radio resolution sample 

^New redshift from Perlman, Schachter & Stocke (1998) 

^Redshift of 0.135 reported in Nass et al. 1996 is incorrect (Bade, private communication) 

References to Table 2: 

B94: Bade et al. , 1994 

B95: Brinkmann et al. , 1995 

B97: Brinkmann et al. , 1997 

B98: Bade et al. , 1998 

B99: Beckmann, 1999 

F93: Fleming et al. , 1993 

F98: Fischer et al. , 1998 

G93: Ghisellini et al. , 1993 

H97: Henstock et al. , 1997 

K92: Kollgaard et al. , 1992 

K96: Kock et al. , 1996 

L96: Lamer et al. , 1996 

L97: Laurent-Muehleisen et al. , 1997 

L98: Laurent-Muehleisen et al. , 1998 

M96: Marcha et al. 1996 

N96: Nass et al. , 1996 

P96: Perlman et al. , 1996b 

PG95: Padovani h Giommi, 1995 

R99: Rector, Stocke h Perlman, 1999 

U96: Urry et al. , 1996 

V96: Vermeulen et al. , 1996 

VLA: VLA Calibrator List 



Notes to individual sources in Table 2: 

0007+472: Optical position given in Kock et al. (1996) is 2.1" from radio position; V=18.3mag 
0040+408: Padovani & Giommi (1995) report object as candidate 
0110+418: Br4ooo=32%; 0(Total)=17.3 mag 
0115+253: ID based on 80 A resolution spectrum 
0152+017: Br4ooo=29%; O (Total) =16.1 mag 

0153+712: Marcha et al. (1996) report P=3%, Br4ooo=33% and WA([Ha,[NII])=17 A; E(Total)=15.5 mag 
0202+088: ID based on 80 A resolution spectrum 

0214+517: Marcha et al. (1996) report P=2% and Br4ooo=15% and Wa<3 A; E(Total)=16.5 mag 

0227+020: ID based on 80 A resolution spectrum 

0314+247: Wa([OII]A3727)=14 A; Wa([OIII]A5007)=12 A; Br4ooo=30%; 0(Total)=17.3 mag 

0450+450: The X-ray emission in the vicinity of this source is highly complex and it is unclear how much (if 

any) of the X-ray emission is actually associated with the BLLac; V=18.5mag 

0650+250: V=15.3mag 

0654+427: Marcha et al. (1996) report P=2-3%, Br4ooo=17% and Wa(Hq;,[NII])=8 AA;0(Total)=17.6 mag 
0656+426: Br4ooo=36%; Star-forming Elliptical?; 0(Total)=15.1 mag 
0710+591: Br4ooo=22%; 0(Total)=17.4 mag 
0916+526: ID based on 80 A resolution spectrum 
0929+502: ID based on 80 A resolution spectrum 

0930+350: Henstock et al. (1997) report source is weak, but no emission lines observed 
1000+225A: Optical counterpart between two radio sources; Upper limit reported to core radio flux 
1012+424: Optical spectrum contains unidentified emission lines with Wa=10A & 4A 
1047+546: Padovani & Giommi (1995) report object as candidate 

1053+494: Rector, Stocke & Perlman (1999) report Wa(Hq!)=2.5 A&; Br4ooo=32%; 0(Total)=16.4 mag 
1136+701: V=15.0mag 

1341+399: ID based on 80 A resolution spectrum 
1427+541: Br4ooo=39%; O (Total) =15. 8 mag 
1428+426: Br4ooo=l%; O (Total) =16. 9 mag 
1454+514: Flat Fielding error at ~8380 AA? 

1516+293: Br4ooo=32%; Star-forming Elliptical?; 0(Total)=17.1 mag 
1532+302: Br4ooo=29%; 0(Total)=14.6 mag 
1534+372: WA([Ha,[NII]])=4 A; Br4ooo not measured 
1555+111: V=15.0mag 
1602+308: E=17.5mag 



1621+377: Strongest line in spectrum in Kock et al. (1996) is Wa=9 A 

1725+118: V=15.8mag 

1743+195: V=16.8mag 

1745+398: Br4ooo=17%; 0(Total)=17.8 mag 

1750+470: Br4ooo=29%; O (Total) =18.1 mag 

1806+698: V=14.4mag 

1813+317: V=17.4mag 

1841+591: Br4ooo=17%; O (Total) =19. 2 mag 

2039+523: Padovani &: Giommi (1995) report object as candidate 
2145+073: Br4ooo=10%; O (Total) =18. 8 mag 

2250+384: Br4ooo=7%; Star-forming EUptical?; E(Total)=14.8 mag 
2319+161: E=17.0mag 

2322+346: Br4ooo=33%; Star-forming Elliptical?; Uncorrected O(Total) 



mag 



o 



